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Abstract 

In this article we quantify the possible effects of the breakdown of Migdal's 
theorem on the electron-phonon coupling constant A, on the critical temper- 
ature T c and on the superconducting gap A by examining different kinds of 
superconducting materials either with low and high critical temperature. We 
use the theoretical approach developed by Grimaldi, Pietronero and Strassler 
[PRB 52, 10516 & 10530, 1995] on experimental data taken both from liter- 
ature and from our recent break-junction tunneling experiments in P^S^Ca 
Cu^Os+x- The results show that a large violation of the Migdal's theorem (as 
in Bi2Sr2Ca C^Os+z) yields to a large increase of A and, in a first approxi- 
mation, to a large but different increase of T c and A. The same theory gives 

no modifications when applied to low-T c conventional superconductors. 
PACS numbers: 74.50.+ r; 74.70.Vy 
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The high-T c superconductors (HTS), included the fullerenes, are all characterized by a 
very low Fermi energy, here Ep, that is comparable with the Debye phonon frequency: this 
involves the breakdown of the Migdal's theorem for the electron-phonon (e-p) interaction 
|jlj], and requires a generalization of the many-body theory of superconductivity. 

Recently this generalization has been made [@-§| by using an extremely simplified model, 
isotropous, tridimensional, mono-band and with the spectral function for the e-p interaction 
represented by an appropriate Dirac function f7|: a 2 F (uj) = 0.5- \-u -5 (u — u>o) where u = 
u\ og = exp (| / +oc duj - F ^ In uj , A = 2 J +o ° du a F ^ and a 2 F (u) ex is the e-p spectral 
function experimentally determined by quasiparticle tunneling measurements. In addition, a 
perturbative scheme for what concerns the parameter (Xuto/Ep) and a not perturbative one 
for A has been used. The main result is that the vertex function shows a complex behaviour 
as function of the phonon wavevector q and frequency u and, in particular, it assumes 
positive values in correspondence of small q. The region where q is low can be favorite by 
the presence of electronic correlations and by effects due to the electronic density of states 
||. This fact makes it possible a strong increase of the T c value, in comparison with the 
results of conventional BCS theory. This approximate vertex-corrected Migdal-Eliashberg 
theory, required by the breakdown of Migdal's theorem and performed up to now only in 
s wave, at T = T c and with li* = 0, gives rise to very interesting phenomena but, in order 
to examine them in detail, it is necessary to solve numerically the Eliashberg equations. 
This can be done by introducing a cutoff value for the wavevector q c , also useful in order to 
obtain some analytic results such as the extension to this theory of the McMillan's formula 
for the calculation of T c As we already pointed out, the main result is the strong T c 

increase which occurs together with a relatively low q c and almost ordinary values for the 
coupling constant (A ~ 0.5 — 1). A complex behaviour is present for the isotopic effect which 
can be very high in correspondence of some parameters' regions, but can also cancel when 
ujq > Ep. It is also interesting to note that intermediate coupling constant values (A ~ 1) 
can reproduce properties which, in the standard Eliashberg theory, correspond to very high 
values (A ~ 3). The theory produces also effects on the normal-state properties, which can 
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draw away from the usual Fermi liquid behaviour. At T = T c , the new Eliashberg equations, 
that take into account the first-order vertex corrections, are MM: 



A(u n )Z(u n ) = ttT c J2 



A(u„ 



{u n - u m ) + tol \uJ r , 



Xa^u, Wmi 9c, A) — arctan 

7T 



(1) 



[1-Zl 



E 



X z (u n , w m , g c , m, A) - arctan 



7T 



(2) 



where A(u n ) is the gap function, Z(u n ) is the renormalization function, u n = (2n — ljirksT 
with n = 0, ±1, ±2, ... are the Matsubara frequencies and m = Uo/Ep is a parameter that 
represents a sort of indicator for the breakdown of the Migdal's theorem. The functions Aa 
and Xz are defined as: 



X A (uj n , uj m , q c , m, A) = A [1 + XP c (u) n , u m , q c , m)\ 

+2X 2 P v (uj n ,uj m ,q c ,m) 



(3) 



X z (u n , u m , q c , m, A) = A [1 + XP v (u n , u m , q c , m)\ 



(4) 



where P c (u n , u m , q c , m) and P v (u n , u m , q c , m) are cumbersome functions defined in the orig- 
inal papers pM, while A is the bare e-p coupling factor not renormalized by the violation 
of the Migdal's theorem. 

Following the usual procedure , it is possible to simplify Eqs. 3 and 4 by calculating 

them at u n = and u m = ujq and then to obtain an approximate expression for T c : 

1.13a;o 



T c (q c ,m,X) 



k By /e(l + m 



exp 



m 



2(l + m) 



exp 



l + X z {q c ,m,X)/ (1 + m) 



(5) 



XA(q c ,m, A) 

In the limit m — > 0, Xz = A a and this formula coincides with the exact result obtained by 



Combescot ]10[ in the weak coupling regime. 



Our original approach consists in numerically solving the system made by Eqs. 4 and 5, 
with T c (q c , m, A) = T^ x , and \z{Qc, m, A) = X ex , where T c ex and X ex are the experimental val- 
ues for the critical temperature and for the e-p coupling constant. In this way we obtain ap- 
proximate values for q c and the bare A value not renormalized by vertex corrections. We use 
the experimental e-p spectral function a 2 F (uj) ex , determined from tunneling experiments, 
to calculate X ex , u Q and the opportune Dirac's function a 2 F (u) = 0.5 • \ ex ■ u ■ 5{uj — ojq) 
that has the property to give the critical temperature closest to the experimental one in 
low-T c superconductors. Of course, the frequency ujq also permits us to calculate the ratio 
m which, as a consequence, is not a free parameter of the model. 

We applied this procedure to four different superconducting materials both low-T c and 
high-T c : Lanthanum, Bismuth, Bai-^K^BiOa and Bi 2 Sr2CaCu208 +x . The first two are 



conventional low-T c superconductors ITTJ characterized by a gap of the order of 0.9-1.3 meV 
and T c pa 5 — 6 K. The main difference among them is the value of the e-p coupling constant 
which is in the intermediate coupling regime for La (A = 0.98) and in the very strong one 
for Bi (A = 2.46). Ba1_.cKj.BiO3 (BKBO) is a well known ceramic superconductor with an 
intermediate T c ~ 24.5 K and a cubic crystallographic structure. Experiments have shown 
that the a 2 F(u> ) can be determined from phonon anomalies in the quasiparticle tunneling 
conductance at eV > A and an intermediate A = 1.23 has been determined [Q. In the past 
years we have largely studied the high-T c superconductor Bi 2 Sr2CaCu 2 08 +x . (BSCCO) for 
the determination of the e-p spectral function ||13|| . Very recently we obtained reproducible 
tunneling results on optimally doped BSCCO single crystals with T c = 93 K fl~4]| . The 
Eliashberg function a 2 F(u ) was determined from break-junction tunneling data and a good 
agreement was obtained at eV > A between the experimental density of states and the 
theoretical one calculated by a direct solution of the Eliashberg equations in presence of 
an energy- dependent normal density of states. The coupling constant was consistent with a 
very strong coupling regime (A = 3.34). Similar results have been recently obtained by other 



groups [I_|. The superconducting properties of the four materials are summarized in Table 
1. It is important to notice that we selected only materials characterized by small or very 



small values of the Coulomb pseudopotential /i* ~ (see Table 1), in order to be consistent 
with the approximate vertex-corrected theory that has been written for /x* = 0. 

Table 2 shows the logarithmic phonon energy ujq calculated by the already mentioned 
equation, the maximum phonon energy of the spectrum u max , the Fermi energy Ep and the 
ratio m = u /E F for the selected superconductors. Actually, in Bi, due to its amorphous 
structure and the consequent particular shape of the spectral function at low energy (i.e. 
for u — > a 2 F (u) ^ buo 2 with b = const) ||11|| , we have used uq = 2A/X where A is the area 
of the spectral function. This ujq value, unlike d;i og , gives the correct critical temperature in 
Bi. Due to the rather small value of u>o (few meV) and the large value of Ep (some eV), the 
breakdown of Migdal's theorem is practically absent in La and Bi (m ~ 10 -3 or less). Quite 
different situation is present in BKBO and BSCCO where u is of the order of 10-20 meV 
while Ep of the order of 100 meV. As a consequence, the violation of Migdal's theorem is 
moderate in BKBO (m ~ 0.14) and strong in BSCCO (m ~ 0.26). 

Figures 1 and 2 show the e-p spectral functions a 2 F{uj) of the selected superconducting 
materials taken from literature [|T1|] or from recent tunneling experiments JTSynj. The tiny 
vertical line at abscissa u Q in every graph indicates the position of the Dirac's function which 
plays the role of ot 2 F(uj) in the approximate vertex-corrected theory. 

As already mentioned, by numerically solving the system of Eqs. 4 and 5 with 

X ex = Xz(q c ,X) 
T™ = T c (g c , A) 

we can determine the cutoff value of the wavevector q c and the bare e-p coupling constant 
A. 

It is important to remark that: 

1) For the comparison with experimental values X ex , we used Xz{m, q c , A) and not 
Aa(^, Qc, A) because the renormalization effect on the coupling constant is present in the 
normal state too; 

2) The bare A values, determined by the above mentioned procedure, are not the results 
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of an exact numerical solution of the system of Eqs. 4 and 5 because such a solution 
generally does not exist. As a matter of fact they are obtained by the average of the two 
values Ai and A2 that fulfil the following conditions: T c (q c , X 2 ) = T^ x , Xz(q c , Ai) = X ex 
and AA = |Ai — A2I is minimum (AXbkbo = 0.08 and AXbscco = 0.3). In practice we 
determine the smallest region in the space of the parameters q c and A that is consistent with 
an approximate numerical solution of the system. In the first three columns of Table 3 the 
results of this procedure are shown for the selected superconductors. 

It is clear that for La and Bi, where the breakdown of Migdal's theorem is absent, we 
obtain A = Xz = Aa = X ex and all the values of q c between 0.2-0.3 and 1 give an approximate 
solution of the system. On the contrary, in BKBO the bare A is reduced of the order of 20% 
with respect to the experimental one and the solution is possible only for an intermediate q c = 
0.6. The largest renormalization is present in BSCCO where the e-p coupling constant the 
superconductor would have had in absence of effects due to the Migdal's theorem violation 
is A = 1.85, about 45% less than the experimental value. In this case the solution is possible 
only for a very small value of the wavevector q c = 0.15. The fourth column of Table 3 reports 
the values of the product m ■ X which is the original parameter of the Migdal's expansion 
and quantifies the amount of deviation from the standard Migdal-Eliashberg theory. 

After having calculated the bare A, we determine a proper scaling factor p that, when 
applied to the a 2 F(u) ex , permits to obtain a scaled e-p spectral function which gives a cou- 
pling constant just equal to A. The values of p are shown in the fifth column of Table 3. This 
scaled a 2 F(uj) is then inserted in a program for the direct solution of the standard Eliashberg 
equations at T = |14| and with p* = to be consistent with the vertex-corrected theory. 



By using this approach we determine a first approximation for the values of the energy gap 
Ax and the critical temperature T c \ the superconductor would have had in absence of the 
breakdown of Migdal's theorem. These values together with the ratio 2Aa//cbT c a are shown 
in the last three columns of Table 3. Of course, due to the practical absence of lowest-order 
vertex corrections, Aa and T c \ of La and Bi are exactly the same as the experimental values. 
In BKBO and BSCCO the effects of these corrections on A and T c are very relevant. Aa 



and T c \ of BKBO are about 33% and 28% lower than the corresponding experimental values 
(see Table 3) giving a ratio 2A^//cb^cA = 4 that is consistent with an intermediate coupling 
regime very similar to the La case. In BSCCO the reductions of Aa and T c \ with respect 
to the experimental values are about 48% and 39%, respectively, giving a renormalized gap 
of 12 meV and a renormalized critical temperature of 56 K. As a consequence of the larger 
renormalization of A in comparison with T c , the ratio 2A\/ ksT c \ reduces to 4.98, a value 
very close to the ratio of Bi (see Table 3). We remark that some conventional low-T c very 
strong-coupling superconductors, like Pbo.sBio.5, have A and 2A/IcbT c values as large as 3 
and 5.12, respectively [TjJ. 

Finally, in Figure 3, we show the dependence of the critical temperature of BSCCO 
on q c and m as determined from Eqs. 3 to 5. The light gray circle on the contour map 
approximately shows the region of the parameters corresponding to the solution of the 
system presented in Table 3. There are two interesting observations: (i) it is clear from Fig. 
3 and it was already stressed in the original papers [|J that, at the increase of the breakdown 
of Migdal's theorem (m — > 0.5), T c increases only if q c is small. For q c of the order of 0.8-1, 
the increase of m produces a reduction of T c ; (ii) the action of removing the effects due to 
the deviations from the standard Migdal-Eliashberg theory, that we have tentatively done 
in the last part of the present paper, does not mean to perform the limits m — > and q c — > 1 
but m — > while q c remains constant (in this case q c = 0.15). As a matter of fact, T c \ ~ 56 
K determined by the direct solution of the Eliashberg equations is consistent with m ~ 10~ 3 
and q c ~ 0.15 as shown by the hollow circle in Fig.3. This guarantees a self-consistency to 
our approach, since it has been calculated that A ~ 2 is compatible with a ratio T c /u>o ~ 0.2 
JTJI which is the value we obtain by using T c \ for the critical temperature. 

In conclusion, the results of this paper show that the breakdown of Migdal's theorem, 
certainly present at various degrees in high-T c superconductors, yields to a relevant increase 
of the experimental e-p coupling constant X ex in comparison with the bare one. This increase 
appears roughly proportional to the degree of violation of Migdal's theorem expressed by 
m ■ A. By solving in direct way the standard Eliashberg equations, we have shown that, in 
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first approximation, T c and A are also increased by the effects of the approximate first-order 
vertex corrections described in Refs. 2-6. The amount of such an increase is greater for 
A than for T c thus leading to an amplified value of the experimental ratio 2A//cbT c . In 
BSCCO, the renormalized values A = 1.85 and 2A A / k B T cX = 4.98 are quite compatible with 
a conventional strong-coupling electron-phonon origin of superconductivity. 

We deeply acknowledge the useful discussions with O.V. Dolgov and C. Grimaldi. 
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TABLES 



Table 1: 







f^ex 


A ex (meV) 


1?(K) 


2A/k B T c 


La 


0.98 


0.039 


0.89 


5.04 


4.10 


Bi 


2.46 


0.091 


1.30 


6.11 


4.90 


BKBO 


1.23 


0.04 


4.5 


24.5 


4.64 


BSCCO 


3.34 


0.0093 


23 


93 


5.74 



Table 2: 





u (meV) 




E F (meV) 


m 


La 


4.49 


15 


3300 


0.0014 


Bi 


2.87 


14 


9900 


0.0003 


BKBO 


14.57 


62 


103 


0.1415 


BSCCO 


21.91 


90 


84 


0.2608 



Table 3: 





A 


Aa 


Qc 


m ■ A 


P 


A x (meV) 


T cX (K) 


2A x /k B T cX 


La 


0.98 


0.98 


0.3-1 


0.0013 


1.00 


0.89 


5.04 


4.10 


Bi 


2.46 


2.46 


0.2-1 


0.0007 


1.00 


1.30 


6.11 


4.90 


BKBO 


1.01 


1.4 


0.6 


0.1429 


1.22 


3.05 


17.71 


4.00 


BSCCO 


1.85 


7.65 


0.15 


0.4825 


1.81 


12.10 


56.34 


4.98 
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FIGURE CAPTIONS 



Fig. 1 Experimental electron-phonon spectral functions a 2 F{ui) of La (upper graph) 
and Bi (lower graph) taken from Ref. 11. The vertical lines indicate the energy positions 
ujq of the Dirac's functions that play the role of the spectral functions in the approximate 
vertex-corrected Migdal-Eliashberg theory. 

Fig. 2 Experimental electron-phonon spectral functions a 2 F(u) of BKBO (upper graph, 
from Ref. 12) and BSCCO (lower graph, from Ref. 14). Details are as in Fig.l. 

Fig. 3 Contour map of the calculated critical temperature of BSCCO as a function of q c 
and m. For details see the text. 
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